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Abstract 



Fermion propagator is computed in a simple model on an extremely 
I anisotropic lattice ^ ^ 1. Fermion determinant is evaluated up to ^""^ 

order. Cliiral condensate is estimated in mean field approximation. 

1 Introduction 

Oh! 

r-| ! Finite-temperature studies of the fermion contribution in lattice gauge the- 

ories were on the agenda over a number of years and both MC experiments 
and analytical models were used to challenge the problem. Calculations with 
dynamical fermions and small lattice spacing, however, are still highly ex- 
pensive. Hence the approximations which would hopefully capture some of 
the essential features of the physics are worth a try. We attempt to study 
the fermionic part of the action (see e.g. 



n=l 



on an extremely anisotropic lattice ^ 1) in a simple model [0, where 
terms in (|1]) proportional to are discarded. Here Uf is the number of 
flavors and ^ = ^ {9,^0) is the anisotropy parameter. The dependance of ^ 
on coupling g and on 'bare' anisotropy parameter C,o = a/ is defined by the 
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condition of independence of physical values on spatialQ a and temporal Qr 
lattice spacings. Dirac operator is defined by 

Dx'x = —^^Uy (x) + ^-^f^J {x') (l + 5°ma^) 5^>^^, 

(2) 

where (x) are link (x, x + v) variables, 7j, are Dirac matrices, and v = 

0. 1, 2, 3. We consider the case of infinitely heavy mirror fermions and choose 
Hamiltonian gauge 

Uo (X, t)^^ = 6aX [1 - (exp {-i(Pa (x)} + 1) 5t',o] , (3) 

where exp {—iipa (x)} = ( Ilt^o"^ (x, t) ] are the eigenvalues of Polyakov 
matrix and Nr is the lattice length in the temporal direction and a, \ = 

1, ...,N are color indices. In particular, for the temporal component of Dirac 
matrix one may write 

D°., = S.... (i^A,„ + 1±2£a;„) . (4) 

with 

Aft (x) = exp {-i^pa (x)} 5t',t~i - (1 + ma^) 5t',f (5) 

For antiperiodic border conditions on fermion fields in temporal direction 
the matrix A^/^ may be diagonalized by discrete Fourier transformation with 
half-integer variables, so 

$f ^ = exp {2m {k + 1/2) {t + 1/2) /N^} (6) 

are the eigenf unctions of the operator Aft- After Fourier transformation we 
get 

Kk' = E *f ^A,,,<|.f = 5,,fc,Afc, (7) 

t,t'=0 

where 

Afc = exp {-i [y?^ - 27r (A; + 1/2)] /N^} - (1 + ma,) (8) 
""^ Further for convenience we put a = 1. 
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are the eigenvalues of Aft- Thereby the inverse matrix is given by 

A^l = iV-i exp {2m {t - t') {k + 1/2) /N,} . (9) 

fc=0 

Common properties of projectors = (^^^)^ and i£2Dl±ie. = g allow to 
write 

. / 1 + 7o 1 - 70 \ 1 + 70 . . X 1 - 70 . . X „x 
/ o ^ + ^ — ^ = o / (^) + o / (^) (10) 



for the arbitrary regular function, in particular the inverse operator (D^^) ^ 
may be written as 

(fl?„)-' = i±^A- + i;^(A-)'. (11) 



2 Ansatz 

To sum over the Fourier variables k in (P) we use a standard trick. Let 
/ (uj) be a rational function of coscu and sincj, which has no poles on lines 
Recj = 27rn. The integral over the whole area of the complex variable uj 
for the periodic function splits into a set of equal integrals over 'elementary' 
bands 2mi < Reu — ^ < 27r(n + l). Due to periodicity, the result of 
integration over any band boundary is 

The loop integral over each band area of uj is equal to zero as well, so ex- 
pressing it as sum of residues one may write 



^/Mtan^^ = 5^re.|/, 



\j)tan^} = 0. (13) 



Taking into account that tan {NrUj/2) has poles at a; = 27r (A; + 12) /Nr, one 
may finally write 



^ Y f (2- (k + /N^l = -i 5^ re. { / (^) tan ^} - $, (14) 



fe=0 ' {tOr} 
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where ujr poles of the / {uj) are located in the band < Re — < 27r. Now 
one may apply the ansatz (|1^ to compute A^j',. In the case considered 



^ ^""^ ~ I c^A'+^O-aV }or -Nr < t - t' 7o ^^^^ 

has poles at ojq = —ipa + imar, so taking into account that / (u) for 
Imtu — i> ±00 we find 

i\ dtt'e''"^" - (Of t + t)e^ m t-t'-i 
(A;l,) = — r iv. , 16 

where T = N~^a~^ and the periodic in tk function Ot-^^t2 is given by (for 

_ j 1 for h > t2 . 
- \ for h<h ■ ^ ^ 

Hence, in our approximation the expression for fermion determinant may be 
found in a closed form from (0) and (p^. 



{t't')u^Xl^ for 0<t-t' <Nr 



3 Fermion determinant 

Since (-D°) ^is known, for the fermion contribution up to the O one can 
write 



' In J d^lj^,d^lj. exp {-Sj. (1/0} = -4'^^^ (1/0 = -S^;"^ (0) - ^S^;"^" (0) r', 



where 

N 



-S'^;"'> (0) = In det (Z^°) = 5^ < 5Z In (cos y^, (x) + cosh ^) - S„ 



X Q = l 



(19) 



with 

-Smirr = Nm/T (20) 

and 

4^^^)" (0) = Sp I {D^y' {D^^y' D"} . (21) 

n=l 
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In case of light {m <^T) fermions we get 

2 \ \ ^^^1 3 



^ \ / / t2ti=0 X n=l 

Re |^7„ (t2, x)^^ <l>„^ (x, n)^^^j^ Ul (ti, x)^^ - (x, n)^^^^^ | . 



:22] 



where 



<^a(x) ^^(x + n) 

tan tan 23 

2 2 ^ ^ 

It is easy to see that Sp^^^" (0) turns into zero in the case of U{1) gauge 
group, if spatial link variables Un are static. 



4 Chiral condensate 



To estimate [ipip) one may, in the spirit of Curie- Weiss method, change all 
ifa (x) for the mean field Ip^ defined by the mean field equation. For example, 
in case of SU(2) gauge group ipi (x) = —ip2 (x) = (p (x) /2 ~ ■^/2. Thereby 
for the chiral condensate we get 



) = T^lnZ = UfT ( ^IndetD 
' dm \ am 



(24) 



Va(x) = </'Q 



or 



/2nf = 1 + ^ 4'^^^" (0) 



m/T 



~'~ 1 + cos {(p/2) 



— - + 0(m2,r'). (25) 



It should be noted that even for infinitely heavy mirror fermion the rem- 
nant term Smirr, defined in (pO]) survives. It is easy to make sure that this 
very term introduces a constant component {2nf) into (ipip), that doesn't 
disappear even at m/T —> 0. The mirror fermion contribution may be totally 
removed "by brute force". If we introduce "the counterterm" = —Smirr 
into the original action (p, but preserve a standard definition for {'^p^|J) , such 
" counterterm" would evidently cancel the undesirable component in the first 



term of 



Therefore, we get {ipip) nj^' 



2 4^^^)" (0) 



m=0 



for m/T 0. 
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5 Discussion 



In the parameter area available for modern computers sea quark contribution 
introduces minor changes in MC data. For example, the comparison of spec- 
troscopy results obtained with dynamical and quenched fermions shows no 
dramatic difference. Apparently at the parameter area of the MC simulation 
the sea quarks simply do not affect the spectroscopy above five to ten per 
cent level 0. Major part of MC data@] helps to trace the sea quark contri- 
bution, rather than locate the area where the sea quarks should obligatory 
be taken into account. 

Suggested model allows to evaluate the fermion contribution analytically 
for arbitrary small values of lattice spacings and infinitely large lattices, where 
MC simulations easily become prohibitively costly. Although such contribu- 
tion can be totally incorporated into the invariant measure and may appear 
almost trivial, it significantly changes the phase structure of the model. Such 
changes, especially in the case of heavy quarks, are very similar to those intro- 
duced by external magnetic field into spin systems. We hope that suggested 
model will help not only to estimate the sea quark contribution at lattices 
with size and spacings unaccessible for MC experiment, but might also assist 
in qualifying the parameter area where the presence of dynamical fermions 
leads to particularly appreciable effects. 
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